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I ' Abstract 

a . 

O ■ In this work we have considered the Taylor series expansion of the dynamic scahng relation 

, , , of the magnetization with respect to small initial magnetization values in order to study 

' the dynamic scaling behaviour of 2- and 3-dimensional Ising models. We have used the 

literature values of the critical exponents and of the new dynamic exponent xq to observe 
the dynamic finite-size scaling behaviour of the time evolution of the magnetization during 
early stages of the Monte Carlo simulation. For 3-dimensional Ising Model we have also 
presented that this method opens the possibility of calculating z and xq separately. Our 
results show good agreement with the literature values. Measurements done on lattices 
with different sizes seem to give very good scaling. 
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1 Introduction 

The determination of the critical temperature and the corresponding critical expo- 
nents for the systems exhibiting second-order phase transitions is the major task 
in many areas of physics. In order to determine the critical parameters of the sys- 
tem there exist two approaches. Traditionally, one uses universality and the scaling 
relations in the equilibrium stages of the statistical mechanical system. In other 
words, one studies the long-time regime of the system. The second approach is to 
study dynamic scaling which seems to exist in the early stages of the quenching 
process in the system exhibiting second-order phase transition. Jansen, Schaub and 
Schmittmann P] have shown that for a dynamic relaxation process, in which a sys- 
tem is evolving according to a dynamics of Model A [2] and is quenched from a very 
high temperature to the critical temperature, a universal dynamic scaling behaviour 
within the short-time regime exists [21 E] . 

Monte Carlo simulations are the most commonly used tools to study statistical me- 
chanical systems in thermal equilibrium around the phase transition point. Creating 
statistically independent configurations is controlled by the autocorrelations in the 
system which depend on the system as well as the simulation algorithm. Around the 
phase transition point, the correlation length grows exponentially and this results in 
the effect known as the critical slowing down. On finite lattices the autocorrelation 
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time diverges as r ~ L^, where L is the hnear size of the lattice and z is the dynamic 
critical exponent. 

Dynamic finite-size scaling relation is a generalized form of the finite-size scaling 
relation valid for the equilibrium stages. For the k^^ moment of the magnetization 
of a system with spatial size L this relation can be written as P] 

M^''\t,e,mo,L) = L'^-^'^'"^ M'^^\tL-\eL^/'' ,moL''°) (1) 

where (5 and v are the well-known critical exponents, t is the simulation time and 
e = (T — Tc) /Tc is the reduced temperature. In the study of the short-time dynamic 
scaling of a system, in addition to the critical exponents, the dynamic critical expo- 
nent z and a new and independent exponent xq appear, where xq is the anomalous 
dimension of the initial magnetization itiq. 

Eq.(0) can be reformulated for the time evolution of the magnetization at early 
stages of the simulation. For very small initial magnetization, time evolution of the 
magnetization undergoes a critical initial increase ^1 E] in the form 

M(t) ~ m^f (2) 
where the exponent 6 is related to the exponent xq by 6* = [xq — l3/i')/z. The 
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power-law behaviour relation (Eq.((21)) of the dynamic finite-size scaling is used to 
estimate the critical exponents 9 and (3 /v. In literature, it has been shown that 
the dynamic behaviour of various spin models obeys the dynamic scaling relations 
(Eqs. (HD and (ED) El 13 El El CDl HD 1121 ; the universality of Xq has also been 
studied and numerical evidence for the existence of the universality of the dynamic 
scaling has been given [3 ■ 

In this work we aim to reconsider dynamic finite-size scaling by approaching the 
problem from a different point of view. Instead of the critical initial increase 
(Eq. Q) of the magnetization used in the literature, we have used the dynamic scal- 
ing relation (Eq.(PQ)) by considering the Taylor series expansion of M(^)(t,e,mo,L) 
with respect to small initial magnetization mo- We have seen that the new ex- 
ponent (xo) exists not only in the initial stages, but also until the effects of the 
initial magnetization of the system diminish. This approach seems to give the crit- 
ical exponents and the dynamic exponents with very high accuracy. In order to 
discuss the proposed method we have considered 2- and 3-dimensional Ising mod- 
els which are well-studied by finite-size scaling techniques for both equilibrium and 
non-equilibruium stages. 

In the next section we present the model, the method and the scaling relation for the 
magnetization. In section 3 numerical results for 2- and 3-dimensional Ising models 
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will be reported. The last section contains the conclusions. 



2 Model and the Method 

In this work we have employed 2- and 3-dimensional Ising models which are described 
by the Hamiltonian 

-l3H^KY,SiSj. (3) 

<i3> 

Here (3 — 1/kT and K — J/kT, where k is the Boltzmann constant, T is the 
temperature and J is the magnetic interaction between the spins. In the Ising 
Model the spin variables take the values Si — :kl. In this work we have performed 
Monte Carlo simulations by using Metropohs algorithm. 

The measurements of the thermodynamic quantities which are obtained using ther- 
malized configurations in Monte Carlo simulations are independent of the initial 
configuration. Here, before the thermalization, the convergence path (to the ther- 
malization) may be different depending on which initial configuration is used. Al- 
though the system reaches the same thermalized configurations in the long-time 
regime, the convergence path is highly sensitive to the initial configuration. In this 
work we propose to consider the dynamics of the thermalization for a system with 
two different initial configurations, in order to obtain the new dynamic index xq. 
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One initial configuration can be considered as the configuration with vanishing ini- 
tial magnetization. The second initial configuration can be taken with arbitrary, but 
small initial magnetization. The dynamic finite-size scaling form of the k^^ moment 
of the magnetization is given by Eq.([Q). At the critical point (e = 0), for small 
initial magnetization (mo) values, the function given in Eq.(^ can be expanded into 
Taylor series around mo = resulting in 



M^''\t,0,mo,L) = L-^^/'^A<(^■)(tL-^0,0)+ (4) 

The first term in this expansion has no mo dependence and it can be calculated 
through Monte Carlo simulations by considering the initial configuration with van- 
ishing initial magnetization. The mo-dependent terms (second and higher order 
terms in the expansion) can be represented by a function Ms{t,0,mo, L) which is 
basically 

M.(t,0,mo,L) = ^^'''^f^J °' "^"^ Lo^o^^oL-"-^^/- + O(mo^) (5) 

and it can be obtained by calculating the difference of M^''\t,0,mo, L) and 
M^'^^t, 0, 0, L). This means it has the form. 



M,(t,0,mo,L) = M(^)(t,0,mo,L) - M^''\t,0,0, L). 
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(6) 



Here M^''\t,0,0, L) corresponds to the first term in the expansion. The terms 
M^'^) (t, 0, mo, L) and M^^^ {t, 0, 0, L) can be measured through Monte Carlo simu- 
lations using two different initial configurations, one with initial magnetization mo 
and the other one with vanishing initial magnetization. 

The function given in Egs.© and (jHI) contains all mo dependence. The thermalized 
configurations are independent of the initial magnetization, hence Ms{t, 0, mo, L) is a 
decaying function. This function is expected to be a function of the autocorrelation 
time and the initial magnetization. In this work, we have chosen an exponentially 
decaying function in the form 

M,(t,0,mo,L) = /(mo,L) e-(*/^(^»'^'\ (7) 

Here t is the simulation time, t{L) is the autocorrelation time and (piL) is an ex- 
ponent which is related to the initial slip. If this form fits the data, two dynamical 
exponents can be calculated separately, /(mo, L) is a mo-dependent function, but 
it is has no autocorrelation time dependence. Hence the new dynamic scaling expo- 
nent Xq can be calculated from the size dependence of this function for small values 
of mo. The autocorrelation time and hence z can be calculated from the decay rate 
of Ms{t,0,mQ, L). The validity of our assumptions has been tested by fitting the 
function given in EqQ to our data and the results are presented in Section 3. 
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For very small initial magnetization values, it is expected that the leading term 
in the expansion (Eq. is dominant and the rest of the terms (represented by 
the function Ms{t, 0, mo, L)) is in the same order of magnitude with the simulation 
errors. Increasing initial magnetization results in a sizeable mo-dependent term in 
Eq. Q. In this range of mo one can expect that Ms{t,0,mQ, L) scales like L^o-l3/i^_ 
Further increase in mo reduces the reliability of such scaling since the contribution of 
higher order terms in the expansion (Eq. (^) becomes effective. In our simulations 
we have used different mo values to investigate the proposed scaling behaviour of 
Ms{t, 0, mo, L). To test the scaling we have used the literature values of P/u and Xq 
for 2- and 3-dimensional Ising models and also calculated the dynamic exponents z 
and Xq for 3-dimensional model through a fit to our Monte Carlo data. 

3 Results and Discussions 

We have studied 2-and 3-dimensional Ising models evolving in time according to 
the dynamics of the Model Aj2]. We have prepared the lattices with varying ini- 
tial magnetization values starting from mo = 0. Totally random initial configura- 
tions are quenched at the corresponding infinite lattice critical temperature. Sim- 
ulations are performed by using Metropolis algorithm on lattices with linear sizes 
L = 512, 640, 768, 896 and L = 32, 64, 96 for 2- and 3-dimensional Ising models re- 
spectively. For error calculations we have binned data points obtained from 4000 to 
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15000 independent runs for 3- dimensional model, and 50000 to 500000 independent 
runs for 2-dimensional model, starting from randomly prepared configurations. For 
3-dimensional model, 3000 iterations are performed for each independent run. For 
2-dimensional model there are two different sets of runs. One set of data is obtained 
after 500 iterations for each independent run and for the other set of data, each 
independent run includes 10000 iterations. 

The dynamic critical exponent z has a little dependence on dimension for local al- 
gorithms. For a given numbers of spins, the saturation in magnetization can be 
observed after relatively small number of iterations in 3-dimensional lattice com- 
pared two 2-dimensional one, since the saturation is observed when the correlation 
length reaches the hnear size of the lattice. In order to observe the behaviour of the 
mo-dependent terms, one needs to study the time development of the system until 
both M('^)(t, 0, mo, ^) and M('^)(t, 0, 0, L) reach approximately the same value. In 
our simulation studies, we have observed that for 3-dimensional lattice, configura- 
tions with both vanishing and non-vanishing initial magnetization reach the same 
value in (approximately) 3000 iterations for the largest lattice (L = 96). For 2- 
dimensional model, since much larger linear sizes are used, the initial configuration 
dependence remains after 30000 iterations and observing the decaying behaviour of 
Ms(i, 0, mo, I/) becomes extremely expensive. For this reason, we have tested our 
ideas for 3-dimensional model. In 2-dimensional lattice, we have performed simu- 
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lations up to 10000 iterations, but even this many iterations is extremely short to 
observe the decaying behaviour of Ms{t,0,mQ, L). Nevertheless, it is sufficient to 
test the scaling behaviour and to obtain the new dynamic critical exponent Xq. 

In Figure 1, the magnetization data are presented for vanishing initial magnetization 
and for ttlq ^ 0.1 in 3-dimensional Ising Model with linear sizes L = 32, 64, 96. The 
magnetization data converge approximately to the same value for both initial con- 
figurations after approximately 500(L = 32), 1500(L = 64), 3000(L = 96) iterations. 
The differences between the data obtained by starting from different initial mag- 
netizations are large initially for large lattices and vanishes after above mentioned 
number of iterations are reached. The final magnetization value decreases as the lat- 
tice size increases, and the magnetization is smaller initially for vanishing rriQ. The 
function Ms{t, 0, mo, L) which represents the difference between M^'^^t, 0, mo, L) and 
M('=)(t,0,0,L) is given in Figure 2, for mo ~ 0.1 and for linear sizes, L = 32, 64, 96 
in 3-dimensional lattice. The function Ms{t,0,mQ, L) is a decreasing function with 
an initial value proportional to the lattice size L. Despite the fact that our main 
aim is not to obtain the dynamical critical exponent {z), since we are interested in 
the functional form of the mo-dependent term and the speed of its decay, we have 
fitted our data to the functional form (r ~ L^) and obtained z from the fit results. 
We have fitted the data to the exponential function given in Eq.Q. This function 
fits the data reasonably well, as it can be seen from Figure 2. The correlation times 
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(r) obtained from the fit are t(L) = 80, 372, 839 for L = 32, 64, 96 respectively. 
The errors are in the order of one. We have fitted the functional form (r ~ L^) to 
our data and we have obtained the value for z as z = 2.04 ± 0.07 for 3-dimensional 
model. Figure 3 shows the fitted function and the data points for three lattice sizes 
used. 

As it is mentioned in Section 2, the function Ms{t,0,mo, L) representing mo- 
dependent terms scales like L^^. Here we have used the literature values of the 
exponent 9 = {xq - P/u)/z, 9 = 0.191(2) [13] and = 0.108(2) llj for 2- 
and 3-dimensional Ising models respectively, to test the scaling. The exponent 
z 9 = Xq — P/u is obtained from these literature values of both 9 and z. For 2- 
dimensional Ising Model, the value of the dynamic critical exponent z for Metropo- 
lis algorithm varies from 2.13 to 2.18 il El UHl CHI I2D1 IHl 122 123 121 ■ For the 
3-dimensional case, the equilibrium case, the value of the dynamic critical exponent 
is 2; = 2.04(3) [25 , while from the dynamic scaling, z ranges from 2.032(11) to 
2.052(17) |TTj. We have used z = 2.18 and z = 2.04 for 2- and 3-dimensional Ising 
models respectively. The scaling of the function Ms{t,0,mo, L) obtained from the 
Monte Carlo data for 3-dimensional model is given in Figure 4. As it can be seen 
from this figure, these terms are scaled reasonably well. We have also fitted the 
data to a function to obtain the best fit and to obtain the value for z9. Our value 
z9 = 0.27±0.03 (for 3-dimensional model) is slightly larger than the value in the 
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literature. 



For 2- dimensional Ising Model, the autocorrelation time (r) is extremely long for 
the lattice sizes we use. Consequently, we didn't aim to see the decaying form of 
Ms{t,0,mo, L). Instead we have tested the scaling behaviour of Ms(t,0,mQ, L). In 
Figure 5, we have presented the data for Ms{t, 0, mo, L) for after 500 iterations and 
Figure 6 shows the scaling of the data when the literature values are used. 

4 Conclusions 

The short-time dynamic behaviour of the magnetization is studied in 2- and 3- 
dimensional Ising models. The time- dependent magnetization profiles with different 
initial magnetization (mo) values exhibit an initial increase during the time scale set 
by the initial magnetization, the dynamic critical exponent [z) and the new dynamic 
exponent (xo). After a long-time {t ~ L^), all profiles gradually join together, since 
in the long-time regime the effects of the initial conditions are ignorable. In this work 
we aimed to study the effects of the initial configurations on magnetization by using 
Monte Carlo simulation. This can be observed by expanding the magnetization 
around mo = 0. This expansion (Eq. ([[J) gives a leading term which is independent 
of mo and the terms containing the information of how the initial configuration is 
prepared. 



11 



In this approach, the dynamic critical behaviour can be studied by using the data 
from the early stages of the Monte Carlo simulation. In order to obtain a relation 
between the size and the initial magnetization dependence of the Monte Carlo data, 
we have considered both vanishing and non-vanishing initial magnetization runs. 
The iteration by iteration differences of these two sets of runs yield the initial mag- 
netization dependence of the system. Using different lattice sizes we have shown 
that the proposed scaling holds. We have also shown that the dynamic finite-size 
scaling is valid not only for the very early stages (within 100 iterations) of the Monte 
Carlo runs, but it is valid until the effects of the the initial magnetization diminish. 

In order to test our proposed approach, we have employed the literature values of 
xq and we have observed a good scaling behaviour for mo-dependent terms 
for the lattice sizes we used. Such a good scaling gives hopes for the applicability of 
this approach to a wide range of statistical mechanical systems for calculations of 
both dynamic and equilibrium critical exponents. 
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Figure Captions 

Figure 1. Magnetization data for 3-dimensional Ising Model with linear lattice sizes 
L = 32, 64, 96, with vanishing and nonzero initial magnetization (mo ~ 0.1) for each 
size, as a function of simulation time t. The curves for vanishing initial magnetization 
start from the origin. 

Figure 2. The function Ms{t, 0, mo, L) corresponding to mo-dependent terms in the 
expansion as a function of simulation time t, for 3-dimensional Ising Model. The 
plot of the function fitted to the form given in Eq. ((Tj) is also shown. 

Figure 3. The autocorrelation time r(L) as a function of the linear lattice size 
L = 32, 64, 96 in 3-dimensional Ising Model. The function fitted to the form t 
is also shown. 

Figure 4. Scaled form of Ms(t, 0, mo, obtained from the simulation data (given 
in Figure 2), for 3-dimensional Ising Model. 

Figure 5. Ms{t, 0, mg, L) obtained from the simulation data as a function of simula- 
tion time t for 2-dimensional Ising Model (mo ~ 0.07). 
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Figure 6. Scaled form of the data given in Figure 5 for function Ms{t, 0, mo, L) for 
2-dimensional Ising Model. 
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